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Mh' Abstract. We present an analysis of the mass and entropy profiles of three poor galaxy clusters (A1991, A2717 and MKW9) 



observed with XMM-Newton. The clusters have very similar temperatures (kr = 2.65, 2.53 and 2.58 keV), and similar redshifts 
(0.04 < z < 0.06). We trace the surface brightness, temperature, entropy and integrated mass profiles with excellent precision 

Q\ I up to ~ 500 h-^ kpc (A1991 and A2717) and ~ 350 h'^ kpc (MKW9). This corresponds to 0.5(0.35) ^,00, where ^200 is the 

radius corresponding to a density contrast of 200 with respect to the critical density at the cluster redshifts. None of the surface 

Cn ' brightness profiles is well fitted with a single j8-model. Double isothermal j8-models provide reasonable fits, and in all cases the 

^ I value of the external /} parameter is consistent with the value found for richer clusters. The temperature profiles have central 

dips but are approximately flat at the exterior, up to the detection limit. The integrated mass profiles are very similar in physical 
units and are reasonably well fitted with the NFW mass model with concentration parameters in the range C200 = 4-6 and 



M200 = 1.2-1.6x10'"* h^Q Mq. A King model is inconsistent with these mass data. The entropy profiles are very similar at large 



\o 
m 

j^— ^ ■ scale, but there is some scatter in the very central region (r < 50 kpc). However, none of the clusters has an isentropic core. 

^b ' We then discuss the structural and scaling properties of cluster mass and entropy profiles, including similar quality 

Q ! XMM-Newton data on the slightly cooler cluster A1983 (kr = 2.2 keV), and on the massive cluster A1413 (kr = 6.5 keV). 

We find that the mass profiles scaled in units of M200 and r2oo nearly coincide, with < 20 per cent dispersion in the radial range 
[0.05 - 0.5] r2oo, where we could compare the profiles without excessive extrapolation. We provide a quantitative test of mass 
profile shapes by combining the concentration parameters of these poor clusters with other values of similar precision from 



Q ' the literature, and comparing with the C200-W200 relation derived from numerical simulations for a ACDM cosmology. The 

data are fully consistent with the predictions, taking into account the measurement errors and expected intrinsic scatter, in the 
C/3 ■ mass range M200 = [1-2 x 10'"* - 1.9 x 10'''] hj^Mg. This excellent agreement with theoretical predictions - a quasi universal 

^ ' cusped mass profile with concentration parameters as expected - shows that the physics of the dark matter collapse is basically 

understood. Scaling the entropy profiles using the self-similar relation 5 oc r, we find a typical scatter of ~ 30 per cent in 
scaled entropy in the radial range [0.05 - 0.5] r2oo- The dispersion is reduced (~ 22 per cent) if we use the empirical relation 
S oc ^"65. The scatter is nearly constant with radius, indicating a genuine similarity in entropy profile shape. The averaged 
scaled profile is well fitted by a power law for 0.05 < r/r2m < 0.5, with a slope slightly lower than expected from pure shock 
heating (a - 0.94 ± 0.14), and a normalisation at 0.1 rioo consistent with previous ROSAT/ASCA studies. These precise XMM 
observations confirm that the entropy profiles of clusters are self-similar down to low mass (kr ~ 2 keV), but that the entropy 
temperature relation is shallower than in the purely gravitational model. This self-similarity of shape is a strong constraint, 
allowing us to rule out simple pre-heating models. The gas history thus probably depends not only on gravitational processes, 
but also on the interplay between cooling and various galaxy feedback mechanisms. 

Key words. Cosmology: observations. Cosmology: dark matter. X-rays: galaxies: clusters, galaxies: clusters: individual: A 
1991, A 2717, MKW 9, Galaxies: clusters: Intergalactic medium 

1 . Introduction Tx) and redshift, z, of the system, such that Q k A(z)T^ (e.g., 

Bryan & Norman 1998 Eke et al. 11998V It has been known 

In galaxy clusters the hot. X-ray emitting gas of the intra- ^ ,. .u . ui u 1 »• a • . • u a 

° -' -' ° ° tor some time that while such relations do exist in observed 

cluster medium (ICM) lies trapped in the potential well of , . ., • . , ,■ . u»i j«- ^t ^ a 

^ ' ^^ ^ clusters, their actual scaling is subtly dmerent from expected, 

the dominant dark matter component. If clusters were formed 

solely through gravitational processes, the properties of dif- ^ith the advent of XMM-Newton and Chandra, observa- 

ferent haloes (e.g.. X-ray luminosity, L^, the gas mass M„,,) ^^""^ ^"^^ numerical simulations are on an equal footing, and so 

would scale with the mass M (or the global X-ray temperature ^^ '^^" ^^^^ *^ modelling of the dark matter collapse and the 

consequent evolution of the ICM as never before. We are be- 

Send offprint requests to: G. W. Pratt, e-mail: gwp@mpe.mpg . de ginning to see some evidence that numerical simulations pre- 
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diet the eon-eet shape for the dark matter distribution from 
~ 0.01 r2oo up to ~ 0.7 r2oo not only in massive systems 
(e.g., Lewis, Buote & Stocke 2003 Pratt & Arnaud 2002 
Pointecouteau et al.,2004i, but also in low mass clusters (Pratt 
& Arnaud I2003> . and that the shape is independent of clus- 
ter temperature (e.g., Arnaud, Pratt & Pointecouteau 2004 1. 
Thus the underlying universality of the dark matter distribu- 
tion, indicated indirectly by the similarity of gas emission 
measure and temperature profiles of relatively hot systems 
(Markevitch et al. .1998 , Neumann & Arnaud 1999 VikhUnin, 
Forman & Jones Q^^^J Irwin & Bregman 2000, De Grandi & 
Molendi 2003 Arnaud, Aghanun & Neumann l2002> is starting 
to be confirmed. 

At the same time, our knowledge of some aspects of the gas 
physics seems wanting, as the departures from the expected 
scaling relations attest. The current consensus is that there is 
some form of non-gravitational process which affects the gas 
and thus modifies the similarity. The non-gravitational pro- 
cesses have been historically divided into either pre-heating, 
where the gas has been heated before being accreted into the 
potential well, by early supernovae and/or AGN activity (e.g. 
Kaiser .1991. Evrard & Henrv lTWTl Valageas & Sil k.l999j . in- 
ternal heating after accretion (e.g. Metzler & Evrard [T994> . and 
cooling (e.g. Pearce et al. 20001. A great deal of theoretical ef- 
fort, including the use of numerical simulations, has been put 
into determining how these processes aff'ect the scaling proper- 
ties of the ICM (e.g. see Borgani et al. 2004 for a review). 

The preferred quantity in which to cast both observations 
and theory is the entropy, (see e.g., Bower [T997> . The entropy 
reflects the accretion history of the gas, but at the same time is 
likely to preserve the imprint of any other physical process. It 
has become traditional to define the 'entropy' as 5 = iiT/nJ , 
which is related to the true thermodynamic entropy via a loga- 
rithm and an additive constant. Pioneering work on cluster en- 



Table 1. Journal of observations. 



tropy profiles can be found in David, Jones & Forman (I1996> . 
In the self-similar framework described above, the en- 
tropy scales with the temperature simply as 5 cc T. Ponman, 
Cannon & Navarro ( I1999> used ROSAT observations to sug- 
gest that the entropy measured at 0.1 r2oo follows the stan- 
dard self-similar scaling at high temperatures but tends to- 
wards a limiting value in the coolest systems. This sparked 
an avalanche of work on non-gravitational heating (but see 
also Bryan 2000, Dave et al. 2002 1. with a particular fo- 
cus on the pre-heating scenario (e.g. Fujita & Takahara 2000 
Bialek et al. "200?, Tozzi & Norman "2001 , Babul et al. 2002, 
Borgani et al. 2002 1: but no definitive picture emerged con- 
cerning the level, actual timescale, or astrophysical source 
of the extra energy injection (Wu, Fabian & Nulsen 1998 
Loewenstein 2000 Kravtsov & Yepes 2000 Bower et al. 2001 



Brighenti & Mathews 120011 Yamada & Fuiita 120011 Nafli & 
Rovchowdhurv l2002t . However, more recent work has shown 
that there is no limiting value to the entropy, and that in fact it 
scales as 5 oc r"^^ (Ponman, Sanderson & Finoguenov l2003> . 
At the same time, better spatially resolved XMM-Newton and 
Chandra observations have shown that the simple prescriptions 
for gas pre-heating, intended to match the limiting value on 
the entropy, were predicting large isentropic cores in groups 
when in fact there are none (Pratt & Arnaud 2003 Sun et 



Cluster 


z 


Rev. 


Mode" 


MOSl 


4p 
MOS2 


pn 


A1991 


0.0586 


584 


FF 


28975 


28878 


19348 


A2717 


0.0498 


558 


FF 


51682 


51647 


44284 


MKW9 


0.0382 


311 


EFF 


29890 


29716 


20720 



Notes: " EPN observation mode: FF = Full Frame, EFF = Extended 
Full Frame; '' Exposure time after flare cleaning 



al. 120031 Rasmussen & Ponman 120041 Khosroshahi, Jones & 
Ponman'2004^. 

The focus has now started to shift from simple pre-heating 
to more advanced models which reflect the energy balance be- 
tween heating and cooling, and the effect of heating on the ac- 
cretion of the gas (e.g. Menci & Cavaliere 2000, Muawong 
et al. |2002; Finoguenov et al. l200ll Xue & Wu l2003l 
Valdarnini 2003 ; Kay, Thomas & TheunsiMn" Voit et al EUU^ 
2003; Voit & Ponman 2003; Kay 2004; Borgani et al. UUU^ . 
These models have need of stringent constraints from obser- 
vations, which XMM-Newton and Chandra are starting to pro- 
vide. 

As shown by Voit et al. J20021 120031 . the two fundamental 
quantities which define the X-ray properties of a relaxed clus- 
ter are the entropy profile of the gas S (r) and the shape of the 
potential well in which it lies. M(r) and S(r) reflect respec- 
tively the physics of the gravitational collapse and thermody- 
namic history of the gas. Low mass clusters are particularly 
interesting targets for the study of the entropy because in these 
systems gravitational and non-gravitational processes affect the 
ICM in roughly equal proportions. In this paper, we present the 
results from high quality XMM-Newton observations of three 
cool clusters, A 1991, A2717 and MKW9, focussing on the gas 
density and temperature profiles (Sect.|3land Sect.|4)'. We de- 
rive the radial distributions of these quantities with good preci- 
sion between ~ 0.01 and 0.5 of the virial radius (r2oo), enabling 
us to calculate the total mass profiles (Sect.|5} and entropy pro- 
files (Sect. |6|l out to the same distance with unprecedented ac- 
curacy. To investigate the scaling and structural properties of 
cluster mass and entropy profiles, we combine these data with 
previously published results for the cooler group A1983 and 
the hot cluster A1413. We will see that the dispersion in the 
scaled mass and entropy profiles is remarkably small. We dis- 
cuss the implications for our understanding of the dark matter 
collapse and the gas specific physics in Sect.0 

Except where otherwise noted, the following cosmology is 
used throughout: Hq - 70 km s"' Mpc"', Q.,„ - 0.3 and Oa - 
0.7, which is denoted ACDMH70. Where we want to examine 
the dependence of scaling on cosmology, we will also use the 
SCDMH50 cosmology: Hq = 50 km s' Mpc', Qn, = 1.0 and 
Oa = 0.0. 

2. Observations and data reduction 

We show the observation details in Tabled For the follow- 
ing analysis we use PATTERNS 1-12 from EMOS event lists, 

' We will discuss the abundance distributions in a companion paper 
(Pratt & Arnaud, in prep). 
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Fig. 1. XMM/DSS overlay images of each cluster. X-ray contours come from the adaptively smoothed, non-background subtracted, [0.3-0.9] 
keV EMOS+EPN image, and are logarithmically spaced. 



and PATTERN from the EPN event list. Each observation 
was cleaned for soft proton flares according to the method de- 



scribed in Pratt & Arnaud (I2002> . using an initial pass in 100s 
bins in a high-energy band (10.-12. keV for EMOS; 12.-14. 
keV for EPN), with a second pass in 10s bins in a broad band 
(0.3-10. keV). To correct for vignetting, the photon-weighting 
method of Arnaud et al. f200ri, as implemented in the SAS 
task evigweight, was applied to each event file. 

We use the dedicated blank-sky event lists accumulated by 
Read & Ponman ("2003 ) as background files, applying the same 
PATTERN selection, flare rejection criteria and vignetting cor- 
rection as described above. Background subtraction of output 
products was undertaken as described in Arnaud et al. (I2002i . 

3. Gas density distributions 

3.1. Morphologies 

The X-ray morphology can give interesting qualitative (and 
quantitative, see e.g., Buote & Tsai I1996> insights into the 
dynamical status of a given cluster. The positions of bright 
serendipitous point sources were obtained from the combined 
EPIC source list in the pipeline products. We then used the 
CIAO utility dmfilth to excise these sources from the im- 
ages and refill the holes with Poisson noise, and then we 
used csmooth to adaptively smooth the images. In Figure ^ 
we show an XMM/DSS overlay of each cluster, where the 
contours come from the adaptively smoothed [0.3-0.9] keV 
EMOS -hEPN image. 

A 1991 and A2717 exhibit symmetric X-ray isophotes, 
suggesting that they are relatively relaxed. The X-ray emis- 
sion of A2717 is centred on its central D galaxy [PL95] 
ACO 2717 BCG, the Brightest Cluster Galaxy as defined by 
Postman & Laurer (■1995J . which has coordinates 00''03'" 12':89, 
-35°56'12'.'3. Note that the SIMBAD coordinates for this clus- 
ter are cuiTently offset by 6'. The X-ray emission of A 1991 
is slightly offset from the optical position of the central D 
galaxy [WCB96] ACO 1991 A (aka NGC 5778)^, coordinates 
14''54"'31^54, -Hl8°38'31'.'l. 



The X-ray emission is of MKW9 is centred on the central 
galaxy UGC 9886, coordinates 15^32"'32n6, -(-04°40'51'.'0, but 
in contrast to the other clusters, the X-ray morphology is not 
symmetric at all. This is especially evident at large scale, where 
the isophotes have a pronounced ellipticity in the NE-SW di- 
rection. At small scale the isophotes are preferentially orien- 
tated in the E-W direction. 



3.2. Surface brightness profiles 

In poor clusters, the emissivity of the hot gas. A, depends sen- 
sitively on the abundance and temperature. As will be seen be- 
low, radial gradients in these quantities exist in all three clus- 
ters, and so in deriving gas density profiles from the surface 
brightness distribution, such variations should be taken into ac- 
count. 

For each cluster, after masking of point sources, the sur- 
face brightness profile was extracted from each event list in the 
[0.3 - 3.0] keV band, with the [0.9-1.2] keV band excluded^ 
Events were binned directly from the event lists in circular an- 
nuli centred at the X-ray emission peak. Background subtrac- 
tion was undertaken separately for each camera; the resulting 
EMOS and EPN data were coadded and then the total profile 
binned such that aS /N ratio of at least 3cr was reached. 

The surface brightness profile was then corrected for emis- 
sivity variations as follows. The temperature and abundance 
values were interpolated to each radius in the surface brightness 
profile by fitting the temperature and abundance profiles with 
functional forms (respectively: an empirical function described 
in Eq. |2l and a lognorm function). The emissivity at each ra- 
dius, A{ff), was then estimated using an absorbed, redshifted 
MEKAL model convolved with the instrument response. The 
surface brightness profile was then divided by A{6) normalised 
to its value at large radii. 



^ In fact, the raw image seems to show a surface brightness edge 
' 10" from the X-ray centroid. These XMM-Newton data do not have 



sufficient spatial resolution to resolve the feature, but the Chandra ob- 
servation (Sharma et al. 2004 1 does show substructures at that scale. 

^ This minimises the contribution from the FeL blend, the feature 
which is the most sensitive to abundance and temperature variations. 
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Fig. 2. The combined EPIC surface brightness profile of each cluster in the [0.3-3.0] keV band, with the [0.9-1.2] keV band excluded and 
corrected for the dependence of the emissivity, A, on radial temperature and abundance variations. Each profile is background subtracted and 
corrected for vigentting. The solid line in each case is the best-fitting BB model convolved with the XMM-Newton PSF; best-fit parameter 
values are given in Table|2| 



Table 2. Results of the double fj (BB) model analytical fits to the gas 
surface brightness profile, errors are 90 per cent confidence. 



Parameter 


A1991 


A2717 


MKW9 


ne.oCh,^"^ cm-3) 


5.61 X 10- 


1.23 X 10-- 


7.86 X 10-3 


re 


1 '42+0-21 


I/QC+O.IS 

^•^^-0.13 


3 '36+ '-08 

-'■•'O-0.66 


P 


0.65!°- 


0.63!°:°? 


7+50 

"■'-0.1 


-^cut 


2:08!°;i 


2;o8!°:l? 


3/40+0.60 

■-'■^"-0.30 


'"c,m 


0' 16+0.01 


0' 3 1+004 


0/35+0.05 
^■■-'■'-0.03 


;rVd.o.f 


197.0/160 


339.3/206 


181.3/133 



3.3. Gas density profile modeiiing 

The corrected surface brightness profile of each cluster (shown 
in Fig. |2} is directly proportional to the emission measure 
profile, EM{r), and can thus be fitted with various paramet- 
ric models for the gas density profile, ne{r). These models 
were convolved with the XMM-Newton PSF (Ghizzardiet al. 
[2001 2002 1 and binned into the same bins as the profile under 
consideration. 

In none of these clusters did we find that a standard y6-model 
was a good description of the entire profile. In all cases the fit 
improves and the reduced ;if^ converges to 1.0 as the central re- 
gions are progressively excluded, which indicates that the outer 
regions are in fact well described with this model. We thus fit- 



ted the surface brightness profiles with the double isothermal 
y6-model described in Pratt & Arnaud J2002I their BB model). 
The fit results are shown in Table |2] Note that the profile ob- 
tained beyond 7?cut in the BB model is the classic y6-model, and 
that in all cases the external/? value is consistent with the mean 



P value for hotter clusters (Neumann & Arnaud[T999| 



4. Temperature distributions 

4.1. Hardness ratio images 

The hardness ratio (HR) of a source is an indirect measure of 
its temperature, and taking advantage of the large bandpass of 
XMM-Newton, we can use the information contained in images 
in different energy bands to get a good idea of the projected 
temperature structure of a given cluster. The method we use 



follows closely that described in Pratt & Arnaud J2003t . ex- 
cept that here we take advantage of the CIAO utility dmf ilth 
to excise and refill with Poisson noise the most obvious point 
sources, and then we use csmooth to adaptively smooth the 
images. 

Source and background images were extracted in the [0.3- 
0.9] keV and [2.5-5.0] keV bands, chosen to avoid complica- 
tions caused by the coupling of temperature and abundance at 
low temperatures due to strong line emission. This is especially 
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Fig. 3. The hardness ratio (HR) image of each cluster is overlaid with the same contour as in Fig.Q These images should be viewed bearing in 
mind that the difference in local background and blank-sky background at low energy has not been corrected for. See the electronic edition of 
the Journal for a colour version of this figure. 



important in these low temperature clusters. We used the total 
EMOS+EPN [2.5-5.0] keV image to define a smoothing tem- 
plate, and subsequently applied this template to each source, 
background and exposure map image in turn. We tried vari- 
ous scales and found that a minimum smoothing scale of 2.5cr 
and a maximum smoothing scale of 4cr was a good compro- 
mise between under- and over-smoothing. The backgrounds 
were normalised based on the effective exposure times and 
subtracted from each image, then the images were divided by 
their exposure maps. The hardness ratio image of each clus- 
ter was then calculated using HR - (image [2.5-5.0] - image 
[0.3-0.9])/(image [2.5-5.0] + image [0.3-0.9]). 

We have not coiTected for the difference between the local 
cluster backgrounds and the blank-sky background at low ener- 
gies, because the effect — being energy-dependent — is rather 
difficult to correct for. This means in particular that i) the ab- 
solute HR values shown in Fig. |3] cannot reliably be converted 
into temperatures, and ii) the decline in HR towards the outer 
regions is an artifact. The latter point is explicitly confirmed in 
Sect. l4.2l below. Despite these caveats, the HR images give an 
excellent idea of the temperature structure, which is the most 
interesting quantity for investigating the dynamical state of a 
given cluster. The HR images of A 1991 and A2717 appear rel- 
atively symmetric, with a clear temperature drop in the centre 
of Al 991. On the other hand, the HR image of MKW9 shows 
sub-structure, suggesting, as does the raw image, that it may 
not be fully relaxed. This is further discussed in Sect. 17. II 
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a colour version of this figure. 



4.2. Radial temperature profiles 
4.2.1 . Annular spectral analysis 

Spectra were extracted in circular annuli centred on the X-ray 
emission peak of each cluster. We optimised the annuli by im- 
posing a 5cr detection, after background subtraction, in the [2.0 
- 5.0] keV band; this was possible for all but the most exter- 



nal annuli, the temperature determinations of which should be 
treated with caution. Regardless of the detection significance, 
we additionally set a lower annulus width of 0.'5, to minimise 
PSF effects (i.e., the annular widths were greater than or equal 
to the diameter enclosing ~ 70 per cent of the energy for the 
on-axis PSF). 



G.W. Pratt & M. Arnaud: Three poor clusters observed with XMM-Newton 



5r 



>^ Projected 

O PSF corrected 

A Deprojected 

D Deprojected PSF correbted 



A199- 



> 

0) 



0) 

ex 




4 
Radius 



10 



Fig. 5. An illustration of the effects of PSF blurring and projection. 
The projected temperature profile of A1991 (stars) is shown here com- 
pared to the PSF-corrected profile (diamonds), the deprojected profile 
(triangles) and the PSF-corrected, deprojected profile (squares). The 
effects play a subtle role in the central regions, but the exterior temper- 
ature extimates are fully compatible with the projected profile. Errors 
are \cr confidence level. The electronic edition of the journal contains 
a colour version of this Figure. 



The EMOS and EPN spectra of each cluster were fitted with 
an absorbed MEKAL model. The absorption was fixed at the 
galactic A^h value in the direction of each cluster'*, as measured 
in the HI survey by Dickey & Lockman ( 1990 1. The free pa- 
rameters of each fit are thus the temperatures and abundances 
(measured relative to Anders & Grevesse 1989 1; the EMOS and 
EPN spectra were fitted simultaneously with the temperature 
and abundance linked. The EPN normalisation was not linked 
to that of the EMOS cameras (it is typically slightly lower). 

We show the resulting projected abundance profiles in 
Fig- S (the temperature profiles are discussed further below). 
There is strong line emission from several elements at the av- 
erage temperature of these clusters, which will be investigated 
in further detail in a future paper. 



4.2.2. Correction for Projection and PSF effects 

The projected temperature profiles show evidence for central 
cool regions of varying strength in all the clusters, the most 
striking example of which is A 1991. It is clear that both PSF 
and projection effects may play a role, and will need to be 
taken into account if we want to recover the best estimate of the 
mass profile. For the following analysis we use the method of 
Pointecouteau et al. (2004i, where the observed annular spec- 
tra, S ^{E), are modelled with a linear combination of absorbed 
isothermal MEKAL models: 



SfiE) = WABS(A^^)^fl,jMEKAL(r;-,Zj). 



(1) 



j=i 



For pure PSF correction, the aij redistribution coefficient is 
the emission measure contribution of ring j to ring /, where the 
coefficients were calculated using the emission measure profile 
from the best-fitting gas density profile model, convolved with 
the XMM-Newton PSF at 1 keV. The redistribution coefficients 
for pure projection effects are calculated as the emission mea- 
sure contribution of the shell j to the ring /. PSF and projection 
effects are taken into account using the emission measure con- 
tribution of the shell j to the ring / after convolution with the 
PSF 

For each cluster the annular spectra were fitted simultane- 
ously in XSPEC. We fixed the absorption to the galactic value, 
and we froze the abundance of each MEKAL model to the best- 
fit value found for each projected annulus. For n annuli, the 
free parameters are thus the n temperatures and n normalisa- 
tions, one per annulus, the other normalisations being linked 
according to Eq.^ Further details of the fitting procedure can 
be found in Pratt & Arnaud J2002t and Pointecouteau et al. 
( {2004 ). 

Figure|5lillustrates the separate effects of PSF blurring and 
projection on the the temperature profile of A 1991 (results for 
A2717 and MKW9 are similar). There is a definite effect in the 
central regions. While in both cases the temperature of the cen- 
tral bin is decreased, the PSF coiTection and deprojection have 
an opposite effect on the second bin. Beyond the third annulus 
(~ 2') the corrected profiles are consistent with the projected 
profile. This is expected, the external profile being very flat. 
However, we can see that the correction process amplifies con- 
siderably the sUghtest variation in the profile and the corrected 
profile is noisier. 

The noise increase is particularly marked for the depro- 
jected, PSF corrected profile, where the largest deviations cor- 
respond to projected annuli which have slightly lower or higher 
temperatures than the best fit smooth profile. Similar behaviour 
was also seen in A2717, MKW9 and in the analysis of A478 



by Pointecouteau et al. J2004> . These deviations are very likely 



* If the absorption is left as a free parameter, the absorption values 
for A2717 and MKW9 compare favourably with those measured in 
the HI survey by Dickey & Lockman <1990> . but that for A1991 is 
significantly higher. However, since leaving the absorption free has 
no discemable effect on the measured temperatures and abundances 
of any of these clusters, we have chosen to fix the Nu to the galactic 
value. 
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Fig. 6. Cluster temperature profiles. These composite profiles are 
combined from the inner three annuli of the deprojected, PSF cor- 
rected profile, plus the projected temperature values thereafter (see 
discussion in Sect. l4.2.'2l . Errors are Icr. The solid lines are the best 
fits to each profile with a function of the form given in Eq. |2l fit pa- 
rameter values are given in Table|3| 



an artifact of the correction process, and as such should not 
be considered as physically meaningful. In the central regions, 
the deprojected, PSF corrected temperatures are subtly differ- 
ent from the projected values. 

Since for all the clusters the outer regions of the depro- 
jected, PSF corrected temperature profiles are i) subject to un- 
physical jumps which would lead to mass discontinuities, but 
ii) consistent with the projected profile, for the following analy- 
sis we use a composite temperature profile made up of the inner 
three annuli of the deprojected, PSF corrected profile, plus the 
projected temperature values thereafter. These composite pro- 
files are shown in Fig.|6l 



Table 3. Results of the best fit of the deprojected, PSF-corrected tem- 
perature profiles (see Sect. 14.2.21 with the analytical model given by 
Eq.ia 



Cluster To 

(keV) 



(keV) 



1 



A1991 1.54 

A2717 1.57 
MKW9 1.76 



1.12 
0.88 
0.72 



0:52 5. 
0:52 2.28 

i:oo 5. 



Table 4. The temperature and abundance values found after simulta- 
neous fits of the EMOS and EPN global spectra of each cluster with 
an absorbed MEKAL model. The global spectrum was extracted in 
the radial range 0.1 ^200 < r < 0.3 r2m, where ^200 comes from the 
best-fitting NEW model discussed later in the text. The abundances 
are relative to Anders & Grevesse >1989j . Errors are Icr for one inter- 
esting parameter. 



Cluster 



kT 
(keV) 



Z 

(Ze) 



A1991 

A2717 
MKW9 



2 65+"-"^ 33+'"" 
2.53+°-«^ 0.34+''02 



^-0.05 
T S9+0.15 



^-0.02 
"■■-"-0.06 



4.2.3. Temperature profile modelling 

We also modelled the composite temperature profile of each 
cluster discussed above (corrected for projection and PSF ef- 
fects) with the function described in Allen et al. J2001> . viz., 

T^To + Td(r/r,)V(l+(r/r,m, (2) 

The parameters for each cluster are given in Tabled 

4.3. Virial temperature 

To compare the properties of these cool clusters with those of 
hotter systems (e.g. Sect.|6j it is useful to define a global tem- 
perature, representative of the 'virial' temperature. Faced with 
the considerations i) that the measurements must be easily re- 
producible both for observers and simulators alike, ii) that in 
none of these observations are we able to detect emission much 
beyond ~ 0.4 r2oo and iii) that all clusters appear to have some- 
what cooler gas in the core region, we extracted a spectrum for 
each cluster from all events between 0.1 r2oo < r < 0.3 r2oo- 
The r2()o in each case comes from the best-fitting NFW mass 
model to the mass profiles, discussed in more detail in Sect. l5.2l 

The EMOS and EPN spectra of each cluster were fit- 
ted simultaneously in the [0.3 - 6.0] keV band with an ab- 
sorbed MEKAL model. Fits were identical to those described 
in Sect. 14.31 The results of this global analysis are given in 
TableEl 

The temperature of A 1991 is in excellent agreement with 
the temperature kr ~ 2.7 keV, derived with Chandra beyond 
the cooling core (Sharma et al. "2004'). The XMM tempera- 
ture of A2717 is higher than the temperature kr = 1.6 + 0.3 



keV derived by Liang et al. J1997> . This value was found from 
an isothermal fit to the ROSAT data within 8', actually an emis- 
sion weighted temperature including the cooling core emission, 
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and is thus likely to be biased low (especially since ROSAT is 
particularly sensitive to any cool component). The global tem- 
perature for MKW9 is in good agreement with previous ASCA 
studies by e.g., Finoguenov et al. (2001 kr = 2.92 + 0.43 keV) 
and Sanderson et al. (^003, kr = 2.88+[J^^ keV). It is also in 
agreement with the analysis of the same XMM-Newton obser- 
vation by Kaastra et al. ^2004 J . 



5. Mass profiles 

5.1. Mass profile calculation 

Combining the gas density (Sect. 13. 3> and temperature 
(Sect. I4.2.2> profiles, we can derive a total gravitational mass 
profile under the assumptions of hydrostatic equilibrium and 
spherical symmetry. The mass was calculated at each radius of 
the temperature profile using the adapted version of the Monte 
Carlo method of Neumann & Bohringer ( 1995 1 described in 
Pratt & Arnaud 



J2003> . This takes as input the parametric 
model for the gas density profile and the measured tempera- 
ture profile with errors. A random temperature is calculated at 
each radius of the measured temperature profile, assuming a 
Gaussian distribution with sigma equal to the Icr error, and 
a cublic spline interpolation is used to compute the deriva- 
tive. Only profiles corresponding to a monotonically increasing 
mass gradient are kept: 1000 such profiles were calculated. The 
input temperature profile is the composite profile discussed in 
Sect. l4T2l 

The uncertainty in the modelling of the gas density profile 
was taken into account by calculation of the errors on the den- 
sity gradient at each point. The final uncertainties on each mass 
profile point are then the quadratic addition of these errors with 
the Monte Carlo errors. 



5.2. Mass profile modelling 

The integrated mass profile of each cluster was then fitted us- 
ing the Navarro et al. C1997. hereafter NFW) density distribu- 
tion p(r) oc [(r/rs)(l + r/^s)]"'. This model has two free pa- 
rameters: a normalisation factor and the scaling radius Ts, or 
equivalently the total mass M200 and the concentration param- 
eter C200 - f2oo/rs (see e.g., Suto, Saski & Makino [19981 for 
details). M200 is defined as the mass contained in a sphere of ra- 
dius r20Q, which encompasses a mean density of 200 times the 
critical density at the cluster redshift pdz) = 3/z(z)^Ho/87rG, 
where h^{z) = f^m(l + z)^ + i^A . In numerical simulations, this 
sphere is found to correspond roughly to the virialised part of 
clusters. Results of these fits are shown in Table|5] Fig.0shows 
the cluster mass profiles plotted in physical units, overlaid with 
the best-fitting NFW models. 

We also tried fitting these mass profiles with an isothermal 
sphere model and a Moore et al. ( II999> model: the former was 
not an acceptable fit to any of these clusters, while the latter 
was unconstrained. 
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Fig. 7. Integrated total gravitating mass profiles, shown with Itr er- 
rors. The solid lines are the best-fitting NFW profiles as detailed in 
Table|5| Dotted lines represent extrapolations of the best-fitting NFW 
models. The data for A1983 and A1413 are from Pratt & Arnaud 
1 2003 1 and Pratt & Arnaud <2002> scaled to the ACDMH70 cosmol- 
ogy. See the electronic edition of the Journal for a colour version of 
this figure. 

Table 5. Results from the NFW fits to the mass profiles. 



Parameter 


A1991 


A2717 


MKW9 


ACDMH70 








C200 
r, ( kpc) 


191!!? 


4-2!!!:^ 

261!g 


5-4t? 
186!« 


rjoo ( kpc) 


1105 


1096 


1006 


M200 (10'*Mo) 


1.63 


1.57 


1.20 


1 1 


9.98/9 


15.8/10 


4.0/8 


SCDMH50 








c 

r. (kpc) 


5 6+0-4 
260:26 


4 1 +"3 
358+37 


5 3+"-^ 
255+*' 


^■200 (kpc) 


1466 


1466 


1358 


M200 (10"Mo) 


2.17 


2.12 


1.63 


x'h 


9.98/9 


15.8/10 


4.0/8 



5.3. Cluster dynamical state and NFW fit 

The NFW model is not a very good fit to the A2717 mass pro- 
file (the;^^ ~ 1.6). Jing (2000 1 has shown that non-equilibrium 
haloes produce poor fits to the NFW distribution, and also that 
the concentration parameters are generally lower in these cases. 
A2717 has the lowest c and the worst-fitting NFW profile, and 
so this may well be due to the halo being unrelaxed. An op- 
tical substructure study by Girardi et al. JI997> has suggested 
that there is a significant structure in the central region, slightly 
foregrounded by ~ 600 km s ', and so the core region may 
indeed be dynamically perturbed. However, we do not see any 
obvious structure either in the X-ray/optical overlay or in the 
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Table 6. A comparison of the dark matter velocity dispersions calcu- 
lated from the a"oM-A^200 relation of Evrard & Gioia <2002l Eq. |3}, 
with the optically-derived galaxy velocity dispersions. Values are 
given in km s"'. 



.000 



Cluster (Tom <^i 



A1991 526 
A2717 520 
MKW9 474 



opl 

631 + 1*' 

541+«] 

579^31 



Reference 
Girardi et al. 11^^ 
Girardi et al. 1 19§3 
Beers et al. «1995i 



hardness ratio map. Moreover, A 1991 and MKW9 are well fit- 
ted with the NFW profile and have similar c values, but while 
A 1991 appears relaxed, MKW9 does not. The isophotal twist- 
ing from the centre to the outer regions (Fig. Q, and the con- 
siderable structure in the hardness ratio image (Fig.^}, together 
with the evidence for substructure from optical studies (Beers 
et al l 1995) . all suggest that MKW9 may not be fully relaxed. 
We can use the tight correlation between M200 and the 
velocity dispersion of the dark matter found in the Hubble 
Volume simulations (Evrard & Gioia 2002 1, viz.. 



ctdm = 1075[/2(z) M2oo/(10'5 h'^, M^)]'/' km s"'. 



(3) 



to further probe the link, if any, between the NFW fits to these 
clusters and their dynamical state. The optically-derived galaxy 
velocity dispersion is a measure of the actual dark matter veloc- 
ity dispersion, both collisionless components presumably fol- 
lowing the same dynamics. In Table |6j we compare the ctdm 
values expected from the best-fitting NFW mass models to 
these clusters with the galaxy velocity dispersions. For unre- 
laxed clusters we would not expect good agreement because of 
i) incorrect estimates of the total cluster mass through the HE 
equation and/or ii) intrinsic departure from the ^^200-0" rela- 
tion. We can see that not only are the results in excellent agree- 
ment^, but that the agreement is particularly good for A2717. 
This gives us confidence in the total mass estimates derived 
from the NFW fits to these X-ray derived mass profiles, and 
suggests that none of these clusters can be too far from equilib- 
rium. 

That being said, it is highly ulikely that the three clusters 
are in exactly the same state of relaxation. Judging by the X-ray 
and optical information, MKW9 is probably the least relaxed. 
In this sample, though, there is no obvious link between the 
dynamical state and the goodness of the NFW fit or the value 
derived for the concentration parameter. 

5.4. Scaled mass profiles 

We now turn our attention to the scaled mass profiles of these 
systems, where we scale the radius by raoo and the mass by 
M200, these values coming from the best fitting NFW model to 
each cluster. Comparison of these scaled profiles allows us to 
assess similarity in the shape of the mass profiles. To increase 
the sample and allow a first comparison with massive clusters, 
we also add data from our previously published work on the 

^ In view of the large error bars on the optical velocity, the agree- 
ment unfortunately does not give stringent constraints on the dynami- 
cal sate of MKW9. 
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Fig. 8. Cluster mass profiles scaled in units of ^200 and M2oo- ''200 and 
M200 were calculated from the best fitting NFW model fSect. l5.2l . See 
the electronic edition of the Journal for a colour version of this figure. 

Table 7. Relative dispersion in scaled mass profiles, measured us- 
ing the standard deviation and mean (first and third columns) or the 
biweight estimators of Beers et al. < 19901 . 



Radius 


ACDMH70 


SCDMH50 




cr/m 


■Sbi/Cbi 


cr/m 


•Sbi/Cbi 


Scaled mass: NFW best fit model 


0.18 




0.05 ^200 


0.18 


0.19 


0.19 


0.1 ^200 


0.15 


0.16 


0.15 


0.16 


0.3 ^200 


0.08 


0.09 


0.08 


0.08 


0.5 rim 


0.04 


0.05 


0.04 


0.04 


Scaled mass: interpolated data 


0.17 




0.05 ^200 


0.18 


0.20 


0.18 


0.1 ^200 


0.25 


0.25 


0.24 


0.24 


0.3 rim 


0.06 


0.02 


0.06 


0.04 


0.5 rim 


0.13 


0.13 


0.13 


0.12 



slightly cooler cluster A1983 (Pratt & Arnaud '2003'i and the 
hot, massive cluster A1413 (Pratt & Arnaud 2002 1. The mean 
temperature in the radial range 0.1r2()o - 0.3r2oo is kr = 2.2 
keV and kr = 6.5 keV for A1983 and A1413, respectively. 
The mass profiles of these two clusters^ are plotted in physi- 
cal units in Fig.Q The scaled mass profiles of all five clusters 
(ACDMH70) are shown overlaid with the best fit scaled NFW 
models in Figure|8] 

It is clear from Fig. 0that the mass profiles of A1991, 
A27 17 and MKW9 are genuinely similar even in physical units. 
When we add the data from A1983 and A1413 and scale the 
profiles, the resemblance is remarkable. To further quantify 
this similarity, we have estimated the relative dispersion in the 
scaled mass profiles at 0.05, 0.1, 0.3 and 0.5 r2oo- The scaled 
masses at these radii were estimated using either i) the best- 



The published data are scaled to the ACDMH70 cosmology. 
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Fig. 9. Cluster entropy (5 = "kTrQ ) profiles, derived from the de- 
projected and PSF corrected temperature profiles and the best-fitting 
analytical model for the gas density. The solid lines represent the en- 
tropy profiles obtained using analytic models for both the gas density 
and temperature distributions. Dotted lines represent extrapolations of 
these models. The data for A1983 and A1413 are from Pratt & Arnaud 
<2003t . See the electronic edition of the Journal for a colour version 
of this figure. 



fitting NFW mass model or ii) a linear inter/extrapolation of 
the measured M{r) in the log-log plane. In the case of A1983 
and MKW9, there is a small amount of extrapolation needed 
to reach 0.5 r2oo- We used both the standard deviation and 
mean and the biweight estimators for location and scale {Cbi 
and S Bi), described in Beers et al. ("1990"l, to estimate the rel- 
ative dispersion. The resulting values are shown in Table 
The scaled mass profiles depend on the cosmological model 
through the function pdz), used in the definition of r2oo and 
M200, and via the angular distance used to convert the angular 
radius to physical units. The cool systems are all at roughly the 
same redshift but A1413 is at a somewhat higher redshift of 
Z = 0.14. We thus also give the dispersions obtained for the 
SCDMH50 cosmology in Table |7] 

As can been seen from Table for the present sample, the 
dispersions are neither sensitive to the cosmology nor to the 
choice of the statistical estimator. The best fitting scaled NFW 
profile of each cluster depends only on the concentration pa- 
rameter. By definition of the scaling the dispersion is null at 
r = '"200, and naturally increases with decreasing scaled radius. 
The observed dispersion is small, less than ~ 15 per cent for 
r > 0.1 r2oo, reflecting the modest dispersion in the concentra- 
tion parameter values (further discussed in Sect. 17. H . Turning 
now to the interpolated data, the dispersion is generally higher 
(e.g. 25 per cent at 0.1r2oo) and its variation with radius is more 
chaotic. This is expected, since in that case statistical eiTors on 
the data contribute to the dispersion. The dispersion obtained 
using the best fitting NFW models is probably more represen- 
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Fig. 10. Scaled entropy profiles. The radius is scaled to r2oo measured 
from the best fitting NFW mass model. The entropy is scaled follow- 
ing the prediction of the standard self-similar model of cluster for- 
mation: 5 oc h(zy*'^T . The scaling is performed using the global 
temperature T, estimated as described in Sect. l4.3l in units of lOkeV. 
See the electronic edition of the Journal for a colour version of this 
figure. 

tative of the intrinsic dispersion in the mass profiles. A more 
rigorous estimate of this intrinsic dispersion would require a 
much larger cluster sample and is beyond the scope of this pa- 
per 



6. Entropy 

6.1. Entropy profiles 

—2/3 

We next determined the gas entropy (S - kTn^. ) profiles 
for the clusters, using the analytic description for the gas den- 
sity profile and the observed temperature profile corrected for 
PSF and projection effects (Sect. I4.2.2> . The resulting pro- 
files, together with the profiles of A1983 and A1413 (Pratt & 
Arnaud 2003 1, are shown plotted in physical units in Fig. |9l 
The profiles are shown with typical errors coiTesponding to the 
eiTor on each temperature bin. We also calculated the entropy 
using the analytical temperature profile model (Eq. |2j. These 
profiles are also plotted in the Figure. 



6.2. Scaled entropy profiles 

In the standard self-similar scenario, clusters form at constant 
density contrast and the gas simply follows the dark matter 
The mean gas and DM density are proportional to the critical 
density of the Universe: n^ '^ Pdm '^ pdz) °^ h^iz). The entropy 
thus scales as 5 oc hizY'^^^T, and the scaled entropy profiles 
{h{z)'^I^T^'^S versus r/r2cn)) of all clusters should coincide. We 
show the entropy profiles with this scaling in Fig. [TO] where 
T is the global temperature estimated as described in Sect. 14.31 
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Fig. 11. Same as Fig[To|but the empirical scaling S cc hizy^^^T^-^^ 
is used instead. This scaling significantly decreases the dispersion 
between the profiles. See the electronic edition of the Journal for a 
colour version of this figure. 

Table 8. Relative dispersion in scaled entropy profiles, measured us- 
ing the standard deviation and mean (first and third columns), or the 
biweight estimators of Beers et al. ( 1990 1. 



Radius 


ACDMH70 


SCDMH50 




cr/m 


Sbi/Cbi 


cr/m 


5'bi/Cbi 


Scaled Entropy: T' 


' scaling 


0.28 




0.05 ^200 


0.30 


0.28 


0.26 


0.1 r2oo 


0.30 


0.29 


0.29 


0.28 


0.3 r2oo 


0.30 


0.29 


0.28 


0.26 


0.5 r2oo 


0.34 


0.33 


0.31 


0.31 


Scaled Entropy: T' 


"■*' scaling 


0.21 




0.05 ^200 


0.22 


0.22 


0.20 


0.1 (-200 


0.22 


0.24 


0.21 


0.21 


0.3 r2oo 


0.20 


0.26 


0.18 


0.20 


0.5 r2oo 


0.24 


0.29 


0.21 


0.19 



in units of lOkeV (hereafter Tiq). Recent results (Ponman et 
al.|2003 1 suggest that entropy scales with temperature such that 
S oc 70 65 rpjjg entropy profiles with this scaling are shown in 
Fig.CH 

As with the mass profiles, we have inter/extrapolated the 
entropy profiles to 0.05, 0.1, 0.3 and 0.5 rioo using the best- 
fitting (analytical) gas density and temperature models for each 
cluster The dispersion in the scaled profiles was once again 
calculated for the two cosmologies using both the standard and 
biweight estimators. The results are shown in Table |S] If the 
standard scaling is used the dispersion is about ~ 30 per cent, 
and is insensitive to the choice of the cosmology or statistical 
estimator There is less scatter in the scaled entropy profiles if 
one uses the empirically-determined scaling of Ponman et al. 
('2003't, S oc T'O-Ss jjjjg jg obvious if one compares Fig.llOland 
Figin The improvement is quantified in Table |8j the typical 



dispersion is ~ 22 per cent, a drop of ~ 40 - 50 per cent from 
that of the standard scaling (the drop is somewhat smaller for 
the biweight estimator). 

The dispersion is remarkably constant with radius (a vari- 
ation of less than 10 per cent), indicating strong similarity in 
the shape of the entropy profile. This is also clear from Fig. II II 
while there is some variation of form in the inner regions, be- 
yond about 0.05 - 0.1 r2oo, the profiles all become approxi- 
mately parallel. This suggests that the form of kr > 2 keV clus- 
ter entropy profiles is not temperature dependent. The shape, 
which is essentially governed by the shape of the density pro- 
file, is not exactly a powerlaw: it flattens slightly with decreas- 
ing radius following the double y6-model. However, within the 
radial range [0.05 - 0.5]r2oo the mean scaled profile is well ap- 
proximated by a powerlaw: 



h(z) 



4/3rr0«5(r) 



470 



0.1^200 



0.94±0.14 



h^o'^ keV cm^ 



(4) 



where the error on the slope is derived taking into account 
the +lcr dispersion around the mean profile. The slope is close 
to, but slightly shallower than, the S oc r'^ behaviour expected 
from analytical modelling of shock heating in spherical col- 
lapse (Tozzi & Norman .2001 J . behaviour which is also seen 



in cosmological simulations (e.g., Borgani et al. l2002t . This is 
further discussed in Sec. 17.21 

The fact that Figures |9l (unsealed profiles) andlllKprofiles 
scaled with the S oc r"-^^ relation and the virial radius) show 
similar scatter is the result of a close coincidence in logarith- 
mic slopes. For r2oo °^ T'-'-^'' (corresponding to M oc r'-^; e.g, 
Finoguenov et al 1200 11 1, and S oc 70 65^ j-j^g j^g 5'_iog r curves 
are translated along a line of slope -0.65/ -0.57 - 1.14, which 
is similar to observed radial variation of S{r) oc 0.94. The 
scaled profile of A1413 (the hottest cluster) is translated a small 
bit lower than those of the cool clusters (since 0.65 > 0.57). 
The effect of scaling with redshift then compensates a small 
amount for this. 

7. Discussion 

7.1. The gravitational collapse of the dark matter 

The observed shape of the mass profile is an important test of 
our understanding of the dark matter collapse. Although the 
exact slope of the profile in the very centre of clusters is still a 
matter of debate, all numerical simulations of structure forma- 
tion predict a universal form with a central cusp (e.g. Navarro 
et al. 120041 and references therein). 

The NFW profile provides the best fit to the three low 
mass clusters studied here, and a King profile (i.e. a profile 
with a core) is rejected. These three poor clusters thus dis- 
play the cusped distribution expected from numerical sim- 
ulations of CDM collapse, in common with A1983 (Pratt 
& Arnaud I2003> and with hotter clusters observed with 
XMM-Newton or Chandra (e.g., David et al. '200 1' Allen et 
al.l2Q0T| Arabadjis et al. 12^021 Pratt & Arnaud 2002. Lewis 
et al. 2003J Pointecouteau et al. 120041 Buote & Lewis l2004t . 
The profiles we have derived are not precise enough for us to 
constain the exact slope in the centre, although we note that the 
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Fig. 12. Variation of the concentration parameter C200 with cluster 
mass M200- The data points are a compilation of the present work 
(A2717, A1991 and MKW9) and data from the literature. We have 
only considered local clusters (z < 0.15) with C200 measurements pre- 
cise to better than ±20 per cent. The data from A1983, A1413, A478, 
A2029 and PKS0745-191 are from Pratt & Arnaud i2003l . Pratt & 
Amaud ^002 1, Pointecouteau et al. <2004> . Lewis et al. <2003> and 
Allen et al. t2001i . respectively. The published PKS0745-191 values 
correspond to a density contrast of 2500, these were converted to C200 
and M200 using the published best fit NFW model. All errors are Icr er- 
rors. Lines: theoretical variation of the mean concentration with mass 
at z = (full line) and z = 0.15 (dashed line) from the numerical sim- 
ulations of Dolag et al. <2004> . The published relation c(M) was trans- 
lated to the virial radius convention used here (see text for details). 
The dotted lines correspond to the mean plus/minus the standard de- 
viation. See the electronic edition of the Journal for a colour version 
of this figure. 



few observations (of massive clusters ) precise enough to con- 
strain this slope favour an NFW type-profile (Lewis et al. l2003l 
Pointecouteau et al. 120041 Buote & Lewis SUDTi. 

There is a remarkable similarity between the mass profiles 
of the cool systems and that of the hot cluster A1413. The dis- 
persion in the scaled profiles is less than 20 per cent in the 
radial range 0.05r2oo - 0.5r2oo, where the various profiles could 
be compared without excessive extrapolation. In other words, 
the shape of the mass profiles from low mass systems up to 
high mass systems is very close to Universal. We now address 
the question of whether this shape, which is defined by the con- 
centration parameter, is quantitatively consistent with the pre- 
dictions. 

Theoretical modeling does not strictly predict a universal 
profile. Lower mass systems should have, on average, higher 
concentration parameters (e.g., Navarro et al. il997i Bullock et 



al. 1200 U . a consequence of their having formed earlier. There 
also appears to be a wide dispersion in the value of c in simu- 
lated haloes of a given mass, probably linked to scatter in their 

To further test the 



formation epoch (e.g. Wechsler et al. I2002t 
theory of the gravitational collapse, we now compare the ob- 
served concentration parameters with the theoretical relation, 
c{M), taking into account the typical dispersion around it. 

Fig. ^] shows the concentration parameter C200 of the 5 
clusters under consideration plotted versus Miao- We have also 
added values from the literature for local clusters iz < 0.15) 
with C2()o measured to a precision of better than ±20 per cent. 
The references used are given in the caption of the figure. 
This sample, although still of modest size, spans a wide range 
of cluster mass: from M2oa - 1.2 x lO''^h:^QM0 to M200 - 
l.SxlO'^h-'Mo. 



Some care is required in comparing to results from numeri- 
cal simulations. The concentration parameter is commonly de- 
fined as the ratio of the 'virial' radius to the scaling radius of 
the NFW model: c - ryir/rs, and the mass as the 'virial' mass 
within rvir. However different conventions for rvh are used in 
the literature: i) ryir - r20() where r2oo is the radius enclosing 



200 times the critical universe density (Navarro et al. I1997t : 
ii) '"vir = '"A, the radius enclosing A times the critical density, 
where A is derived from the spherical collapse model (Eke et 



al. 120011 Bullock et al. .2001. : Wechsler et al. .2002. Zhao et 
al. I2003t : iii) r^^ = ?"20on the radius enclosing 200 times the 



mean universe density (Dolag et al. I2004t . 

These different conventions lead to different concentrations 
and masses for a ACDM Universe since r2oo < '"a < '"aoon 
(see Huffenberger & Seljak 2003 1. When the various results 
are translated to a common virial radius convention, there is 
good agreement on the c(M) relation (see e.g. Zhao et al. 
2003 and Dolag et al. 2004). We thus only consider the re- 



cent work of Dolag et al. J2004> . who found that the mean 
concentration is well fitted by a power law function of the 
mass, (1 + z)c - co(M/Mo)" with a ~ -0.1 (ACDm cos- 
mology, cTg, = 0.9). Note that these authors used the conven- 
tion rvir = '"20on while we use rvir - '"200 ■ We thus converted 
(using the NFW profile) their C2oon and M2oon values to C200 
and M200 values. At z - we found C2oo/c20on ~ 0.6 and 
M20Q/M2Q0ii - 0.73 - 0.7 in the mass range considered here. 
The corresponding C2oo(A^20o) relation is plotted as a full line 
in Fig. [21 Dolag et al. (2004 1 also studied the scatter at a given 
mass: the concentration follows a log-normal distribution with 
a dispersion of 0.2, in good agreement with previous studies 
(Eke et al. 2001 Bullock et al. l2001> . The curves correspond- 
ing to the mean plus/minus the standard deviation are plotted as 
dotted lines in Fig. El The concentration depends on redshift, 
so we have also plotted the C2oo(A^20{)) relation at z = 0.15, 
the highest redshift for our sample (dashed line). The evolution 
effect is negligible as compared to the dispersion. 

It is clear from Fig.^lthat there is an excellent agreement 
between the observed concentration parameters and the theo- 
retical predictions. Taking into account the measurement er- 
rors and expected intrinsic scatter, there is no obvious deviation 



G.W. Pratt & M. Amaud: Three poor clusters observed with XMM-Newton 



13 



from the theoretical c{M) curve^. However, the present data set 
does not allow us to check the expected increase of concentra- 
tion parameter with decreasing mass. The mass dependence is 
small in the mass range under consideration, being of the order 
of the typical dispersion (see Fig. I12> . We could only test the 
expected c(M) dependence with a much larger sample (prefer- 
ably with individual C200 values measured to a precision smaller 
than the scatter). 

7.2. The gas specific physics 

7.2.1 . A universal entropy profile? 

The present study confirms the remarkable self-similarity in the 
shape of the entropy profiles down to low mass (kr ~ 2 keV), 
suggested by the study of Ponman et al. (2003 1 and our pre- 
vious XMM-Newton test-case comparison of the cool cluster 
A1983 and the massive cluster A1413 (Pratt & Arnaud'200?). 
Note that the accuracy of the entropy profiles, which are well 
resolved from 0.01 raoo up to 0.5 raoo, has allowed us to 
perform a direct check of the self-similarity. In contrast, the 
ROSAT/ASCA study of Ponman et al. ( 2003 1 had to rely both on 
extrapolations beyond the detection radius (particularly at low 
mass) and stacking analysis. These XMM-Newton data have 
also allowed us to quantify the shape of the entropy profile, 
which is found to be only slightly shallower than expected in 
the pure shock heating model. 

Our study also confirms that the S-T scaling relation is 
shallower than predicted by the purely gravitational model. 
The dispersion in the scaled profiles is decreased when we use 
the empirical relation S oc t^.gs^ estimated by Ponman et al. 



J2003> at 0.1 r2oo- We note that the normalisation of this re- 
lation is consistent between the ROSAT/ASCA sample and the 
present XMM-Newton analysis. Ponman et al. ( 2003 ) obtained 
S(0.lr2oo) ~ 530 T°f h^^J^^keV cm^ (their Figure 4). The nor- 
malisation is only 13 per cent higher than the value, we derive 

— 1/3 9 

from the mean scaled entropy at that radius (470 h.^^ keV cm ; 
see Eq.|3l. This is well within the typical dispersion. Note that 
the self-similarity of form found here implies that the slope 
of the S-T relation (also studied by Ponman et al. (2003) 
at rsoo) should not depend on scaled radius, at least up to 
0.5r20() (the maximum radius considered here). Its normalisa- 
tion of course does: from Eq.|4] the normalisation at 0.5 r2oo is 
~2100h^g'^^keVcm2. 

An alternative scaling based on preheating models was pro- 
posed by Dos Santos & Dore (2002i, in which the entropy 
should scale according to 5 oc (l-\-T/To), with Tq ~ 2 keV. Our 
present sample does not allow us to distinguish between this 
scaling and that of Ponman et al. (2003). In fact, the difference 
between scaling by 5 oc r°-*^ and scaling by 5 oc (1 + T/Ty) is 
only ~ 10 per cent in the 2 to 10 keV domain, which explains 
why Ponman et al. (2003) found good agreement when their 



' In our published work on A 1983 (Pratt & Amaud l2003> we men- 
tioned that the concentration parameter is lower than expected for a 
cluster of this mass. This is not strictly true, as we in fact compared 
our C200 value with the Ca value of Eke et al. <2001> for a ACDM 
cosmology, and did not take into account the expected theoretical dis- 
persion. 



entropy profiles were scaled by the latter factor (their Fig. 2). 
We further note that the Dos Santos & Dore scaling is based on 
a preheating model where the factor To was fixed to recover the 
limiting entropy value proposed by Lloyd-Davies et al. ('2000>. 
We now know that such models are inconsistent with both the 
observed power law S-T behaviour and the lack of constant 
entropy cores in cool systems. 

Our sample is still of modest size, and the dispersion analy- 
sis on the scaled entropy profiles should not be overinterpreted. 
It is a first attempt to quantify the degree of self-similarity, as 
such better than a simple comparison of the profiles by eye. 
While the analysis has begun strongly to suggest self-similarity 
at low mass, our conclusion on the self-similarity between poor 
and hot clusters still relies on the assumption that A1413 is 
a fair representative of massive clusters. The dispersion study 
does show that the S oc 7"**^^ scaling is better than the standard 
scaling, but does not allow us to gain any further constraints 
on the exact slope of the S-T relation (our present sample be- 
ing too small, particularly with only one cluster at high tem- 
perature). We recall that Ponman et al. (2003) give a slope of 
a = 0.65 + 0.05. 

We feel it should also be pointed out that the degree of sim- 
ilarity emphasised here does not imply a strict scaling of all the 
gas properties. 

Firstly, there is large dispersion in the central region, r < 
0.05r2oo, roughly coiTesponding to the cooling core region. For 
instance the entropy of A1991 and A1983 at 0.01 r2oo differ by 
an order of magnitude. These huge differences, also observed in 



Chandra data by Sun et al. ( I2004t . are likely to reflect the large 
variety of cooling core histories, such as would be expected 
from the complex interplay between cooling and central AGN 
activity. 

Secondly, relatively modest differences in entropy translate 
into much larger differences in other quantities, such as the X- 
ray luminosity. This is because the luminosity is much more 
sensitive to the gas density than the entropy: Lx oc nl while 
S oc Hg ' , so that Lx oc 5 "-^ for clusters of the same tempera- 
ture and internal structure. This has important consequences for 
sources of the scatter in the Lx-T relation, even after the exclu- 
sion of the cooling core region. To further illustrate this point, 
we can compare in greater detail the 3 poor clusters A 1991, 
A1983 and MKW9, which have practically identical tempera- 
tures (see Tab.|3}. Above 0.1 r2oo (~ 2'), safely outside the cool- 
ing core region, the temperature profiles are isothermal (Fig.|6j 
and the difference in the entropy profiles reflects differences in 
the gas density profiles. MKW9 has a higher entropy through- 
out its ICM, at least up to the detection limit, than either A 1991 
or A2717, by a factor ~ 1.4 - 1.7 (FiglHll- Note that this dif- 
ference corresponds to ±(1. - 1.5)cr around the mean. From 
the simple argument, Lx °^ 5"^, we expect the luminosity of 
MKW9 to be considerably lower (by a factor 3-5) than the 
luminosity of A 1991 or A2717. A direct computation of the lu- 
minosity of the three clusters shows that this is indeed the case. 
The luminosities of A2717 and A 1991 within the virial radius, 
excluding the central 0.1r2oo region, differ by less than 20 per 
cent, while the luminosity of MKW9 is 4 times lower. 

In spite of the above caveats, the picture that emerges from 
the present XMM-Newton study, together with previous con- 
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straints on the slope of the S-T relation (Ponman et al. l2003> . 
is that; 

- There is no break of self-similarity in the entropy profiles 
down to kr ~ 2keV and in the radial range 0.05 raoo < r < 
0.5 r2oo 

- In that radius and temperature domain, the entropy typi- 
cally behaves as Sir) oc ^"■''5±°•°5(r/r2oo)°•'''*='°■'^ with a 
normalisation S ~ 500 h^^' keV cm^ at T = 10 keV and 
r - 0.1 r2oo- 

- The best guess Icr dispersion in the scaled entropy profiles 
is ~ 25 per cent. 

7.2.2. Comparison with theoretical expectations 

The current view is that the departures of the gas proper- 
ties from the standard self-similar picture are due to non- 
gravitational processes. However, in principle, these departures 
could also be due to some flaw in the modeling of the gravita- 
tional collapse itself, e.g. the underlying dark matter compo- 
nent does not obey self-similarity. We would like first to em- 
phasize that the present study, which provides the first precise 
estimates of mass profiles for a significant number of low mass 
clusters, demonstrates that this is not the case (see Sect. 17. 1> . 
However, this does not mean that the gravitational heating of 
the gas is perfectly understood. For instance, Valegeas et al. 
1 1200 3^ suggested that shock heating at very large scales might 
contribute significantly to the cluster entropy. 

It is now fairly clear that pure cooling or simple pre-heating 
models fail to explain the observed cluster properties. This 
has already been greatly discussed in the literature (e.g. see 
Ponman et al. I2003> . and we will only briefly comment on 
it. 

Spherical pre-heating models predict both a break in 
the S-T relation and large isentropic cores (e.g Tozzi & 
Norman 2001 Figures 5 and 17) that are simply not seen, as 
the present study confirms. This does not mean that pre-heating 



Thus the current consensus is that some combination of 
cooling and feedback acts to modify the entropy in clusters. 
The present entropy profile results can be compared with the 
predictions of the recent large scale numerical simulations of 



does not play a role. As pointed out by Voit et al. ( 12003> and 
Ponman et al. (2003'l, pre-heating may also aff'ect the genera- 
tion of intracluster entropy by smoothing the accreted matter 
distribution, an effect not taken into account in simple spher- 
ical models. This effect, in combination with cooling, could 
explain cluster entropy properties (Voit & Ponman 2003 1, al- 
though cosmological simulations are required to fully assess 
this explanation. 

Pure cooling models are not actually discrepant with the 
entropy properties outlined above. They can explain the S-T 
relation at O.lraoo (Muawong et al. 2002" Fig. 4; Kay "2004" 
Fig. 1 ; see also Voit et al. 2002) . and do not predict a strong 
break of self- similarity of the entropy profiles (except perhaps 
below kr ~ 1 keV, Dave et al. 120021 Fig. 5), although they may 
slightly overpredict the entropy at very large radii ("Kav 2004 1. 
The main problem with pure cooling models, as shown by 
the authors, is that they tend to suffer from overcooling (e.g., 
Muanwong et al. 2002 Dave et al. 2002 1. which is at odds with 
the observed mass fraction of the stellar component (Balogh et 
al."200ri. Furthermore, it seems plausible that the feedback be 
associated with the cooling (Voit and Brvan 2001t . 



Borgani et al. (2004) and Kay J2004> . Borgani et al.'s simula- 
tions incorporate a physically motivated model of star forma- 
tion and SNII driven galactic winds. The heating was found 
to be insufficient, leading to a steeper S-T relation than ob- 
served. Moreover, although the simulated entropy profiles were 
very similar in shape, their slope is shallower than found in our 
study: behaving roughly as S(r) oc r^-''^ (from Borgani et al. 
|2004, Figure 14). This is linked to strongly decreasing tem- 
perature profiles, which, as noted by the authors, is at odds 
with the observations. Targeted feedback models (Kay et al. 
120031 Kay 2004 1 seem to be more successful in reproducing 
the observed properties of cluster entropy profiles. In this phe- 
nomenological model, a fraction of particles which are about 
to undergo cooling are reheated to a fixed entropy level. When 
this level is fixed to about 1000 keV cm^, this model repro- 
duces the S-T relation at 0. lr2oo (and at 0.5r2oo)- We also note 
the remarkable agreement between these simulations and our 
observations when comparing the entropy profile slope: both 
are slightly shallower than the shock-heating case. For instance 
the entropy profile of a 3keV cluster above 0.05r2oo is essen- 
tially a power law, with an entropy of ~ 100 keV cm^ and 
1000 keV cm^ at 0.05r2oo and 0.5r2oo respectively. This can be 
compared to ~ 120 keV cm^ and 1000 keV cm^ from our re- 
sults (Eq.0}. It remains of course to find a physical motivation 
for this type of feedback. 

8. Conclusion 

We have presented results from new XMM-Newton observa- 
tions of three poor clusters (A1991, A2717 and MKW9) hav- 
ing similar temperatures (kr = 2.65, 2.53 and 2.58 keV), and 
similar redshifts (0.04 < z < 0.06). The paper concentrates on 
the properties of the mass and entropy profiles, which we were 
able to map up to ~ 0.5 r2oo- We then combined these data with 
previously published data on A1983 (kr - 2.2 keV) and on 
the massive cluster A1413 (kr - 6.5 keV), and examined the 
scaling properties of the profiles. The emerging picture from 
these XMM observations is that local clusters do form a self- 
similar population down to low mass (M200 ~ lO''* h^p M© or 
kr ~ 2 keV). 

Our study has provided clear evidence that the dark matter 
profile of local clusters is nearly universal and presents a central 
cusp, as predicted by numerical simulations. The concentration 
parameter of the clusters, and of other massive clusters from the 
literature, were found to be consistent with the C20()-M2oo rela- 
tion derived from numerical simulations for a ACDM comol- 
ogy, in the mass range M200 - [1.2 x 10'"^ - 1.9 x 10'^] h^^MQ 
(taking into account the measurement errors and expected in- 
trinsic scatter). This excellent agreement with theoretical pre- 
dictions shows that the physics of the dark matter collapse is 
basically understood. 

Except in the very centre, the entropy profiles of these clus- 
ters are self-similar in shape, with close to power law behaviour 
in the 0.05 r2oo < r < 0.5 r2oo range. The slope is slightly shal- 
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lower than predicted by shock heating models, 5 (r) oc ^0 94±o.i4 
We have confirmed that the S-T relation is shallower than in 
the purely gravitational model, and has a normalisation consis- 
tent with that found from previous ROSAT/ASCA studies. The 
entropy scaling behaviour is summarized in our Eq. |3 We em- 
phasize that the self-similarity of shape is a strong new con- 
straint, and simple pre-heating models can already be ruled out. 
In addition to the gravitational effect, the gas history probably 
depends on the interplay between cooling and various galaxy 
feedback mechanisms. 

The shape and scaling properties of the mass and entropy 
profiles of clusters are the key observational constraints to be 
considered by theoretical work. These profiles reflect respec- 
tively the physics of the gravitational collapse and thermo- 
dynamic history of the gas. They can now be measured with 
XMM-Newton in both massive and low mass clusters over ex- 
ceptionally wide radial ranges. Larger cluster samples are still 
required to firmly establish the exact slopes and intrinsic scatter 
of the scaling laws, the S-T, M20[)-T and c-M relations, and to 
confirm the self-similarity of form of the profiles over the full 
mass range. 

Acknowledgements. We thank Doris Neumann for her contribution 
to the early part of this project, and she, Etienne Pointecouteau and 
Klaus Dolag for useful discussions. We thank the referee for com- 
ments. GWP acknowledges support from the French Space Agency 
(CNES), and from the European Commission through a Marie Curie 
Intra-European Fellowship under the FP6 programme (contract no. 
MEIF-CT-2003-500915). The present work is based on observations 
obtained with XMM-Newton an ESA science mission with instru- 
ments and contributions directly funded by ESA Member States 
and the USA (NASA). This research has made use of the NASA's 
Astrophysics Data System Abstract Service; the SIMBAD database 
operated at CDS, Strasbourg, France; the High Energy Astrophysics 
Science Archive Research Center Online Service, provided by the 
NASA/Goddard Space Flight Center and the Digitized Sky Surveys 
produced at the Space Telescope Science Institute. 

References 

Allen, S.W., Schmidt, R.W., & Fabian, A.C. 2001, MNRAS, 328, L37 
Anders, E., & Grevesse, N. 1989, GeCoA, 53, 197 
Arabadjis, J.S., Bautz, M.W., & Garmire, G.P 2002, ApJ, 572, 66 
Amaud, M., Neumann, D., Aghanim, N. et al. 2001, A&A, 365, L80 
Amaud, M., Aghanim, N., & Neumann, D. 2002, A&A, 389, 1 
Amaud, M., Majerowicz, S., Lumb, D., et al. 2002,A&A, 390, 27 
Amaud, M., Pratt, G.W., & Pointecouteau, E., 2004, MmSAI, 75, 729 

lastro-ph/0312398i 
Babul, A., Balogh, M.L., Lewis, G.R, & Poole, G.B. MNRAS, 330, 

329 
Balogh, M., Pearce, RR., Bower, R.G., & Kay, S.T. 2001, 
Beers, T.C., Flynn, K., & Gebhardt, K., 1990, Al, 100, 32 
Beers, T.C., Kriessler, J.R., Bird, CM., & Huchra, J.P, 1995, Al, 109, 

874 
Bialek, J.J., Evrard, A.E., & Mohr, I., 2001, Apl, 555, 597 
Borgani, S., Govemato, R, Wadsley, J.,et al. 2002, MNRAS, 336, 409 
Borgani, S., Murante, G., Springel, V. et al. , 2004, MNRAS, 348, 

1078 
Bower, R.G., 1997, MNRAS, 288, 355 
Bower, R.G., Benson, A.J., Lacey, C.G., et al. , 2001, MNRAS, 325, 

497 



Brighenti, R, & Mathews, W.G., 2001, Apl, 553, 103 

Bryan, G.L. 2000, Apl, 544,L1 

Bryan, G.L., & Norman, M.L., 1998, Apl, 495, 80 

Bullock, I.S., Kollatt, T.S., Sigad, Y., et al. 2001, MNRAS, 321, 559 

Buote D., & Tsai, I., 1996, ApJ, 458, 27 

Buote D., & Lewis, A.D., 2004, ApJ, 604, 116 

Dave, R., Katz, N., & Weinberg, D.H. 2002, ApJ, 579, 23 

David, L.P, Jones, C, & Rorman, W., 1996, ApJ, 473, 692 

David, L.P, Nulsen, PE.L, McMamara, B.R., et al. 2001, ApJ, 557, 

546 
De Grandi, S., & Molendi, S., 2002, ApJ, 567, 163 
Dickey, J.M., & Lockman, RJ. 1990, ARA&A, 28, 215 
Dolag, K., Baertelman, M., Perrotta, R, et al. , 2004, A&A, 416, 853 
Dos Santos, S., & Dore, O., 2002, A&A, 383, 450 
Eke, V.R., Navarro, J.R, & Rrenk C.S., 1998, ApJ, 503, 569 
Eke, v., Navarro, J.R, & Steinmetz, M., 2001,ApJ, 554, 114 
Evrard A.E., & Henry J.P 1991, ApJ, 383, 95 
Evrard, A.E., & Gioia, I. 2002, Merging Processes in Galaxy Clusters, 

ed. L. Reretti, I.M. Gioia, G. Giovannini, Astrophysics and Space 

Science Library, Vol. 272, 253 
Rinoguenov, A., Reiprich, T, Bohringer, H., 2001, A&A, 368, 749 
Rinoguenov, A., Amaud, M., & David, L.P 2001, ApJ, 555, 191 
Rinoguenov, A., Borgani, S.,Tomatore, L., & Bohringer, H., 2003, 

A&A, 398, L35 
Rujita, Y., & Takahara, R, 2000, ApJ, 536, 523 
Ghizzardi, S. 2001, EPIC-MCT-TN-OU (XMM-SOC-CAL-TN- 

0022) 
Ghizzardi, S. 2002, EPIC-MCT-TN-012 
Girardi, M., Escalera, E., Radda, D. et al. , 1997, ApJ, 482, 41 
Huffenberger, K.M., & Seljak U., 2003, MNRAS, 340, 1199 
Irwin, J.A., & Bregman, J., 2000, ApJ, 538, 543 
ling, YP 2000, ApJ, 535, 30 
Kaastra, J.S., et al., 2004, A&A, 413, 415 
Kaiser, N. 1991, ApJ, 383, 104 

Kay, S.T, Thomas, PA., Theuns, T, 2003, MNRAS, 343, 608 
Kay, S.T, 2004, MNRAS, 347, L13 
Khosroshahi, H.G., Jones, L.R., Ponman, T.J., 2004, MNRAS, 349, 

1240 
Kravtsov, A.V., Ypes, G., 2000, MNRAS, 318, 227 
Loewenstein, M. 2000, ApJ, 532, 17 

Lewis, A.D., Buote, D.A., & Stocke, J.T., 2003, ApJ, 586, 135 
Liang H., Pierre M., Unewisse A., & Hunstead R.W. , 1997, A&A, 

321,64 
Lloyd-Davies, R., Ponman, T, & Cannon, D.B., 2000, MNRAS, 315, 

689 
Markevitch, M., Rorman, W., Sazazin, C, & Vikhlinin, A., 1998, ApJ, 

503, 77 
Menci, N., & Cavaliere, A., 2000, MNRAS, 2000, 31 1, 50 
Metzler, C.A., & Rvrard, A.R., 1994, ApJ, 437, 564 
Moore, B., Quinn, T, Govemato, R, Stadel, J., & Lake, G., 1999, 

MNRAS, 310, 1147 
Muanwong, O., Thomas, PA., Kay, S.T, & Pearce, RR. 2002, 

MNRAS, 336, 527 
Nath, B., & Roychowdhury, S., 2002, MNRAS, 333, 145 
Navarro, J.R, Rrenk, C.S., & White, S.D.M. 1997, ApJ, 490, 493 

(NRW) 
Navarro, J.R, Hayashi, R., Power, C. et al. , 2004, MNRAS, 349, 1039 
Neumann, D., & Bohringer, H. 1995, A&A, 301, 865 
Neuman, D., & Amaud, M. 1999, ApJ, 542, 35 
Pearce, R.R., Thomas, PA., Couchman, H.M.P, Rdge, A.C, 2000, 

MNRAS, 317, 1029 
Pointecouteau, R., Arnaud, M., Kaastra, J., & De Plaa, J., 2004, A&A, 

in press,_astro-ph/0403596. 



16 G.W. Pratt & M. Amaud: Three poor clusters observed with XMM-Newton 

Ponman, T.J., Cannon, D.B., & Navarro, J.F. 1999, Nature, 397, 135 
Ponman, T.J., Sanderson, A.J.R., & Finoguenov, A. 2003, MNRAS, 

343, 331 
Postman, M., & Lauer, T.R. 1995, ApJ, 440, 28 
Pratt, G.W., & Amaud, M. 2002, A&A, 394, 375 
Pratt, G.W., & Amaud, M. 2003, A&A, 408, 1 
Rasmussen, J., Ponman, T.J., 2004, MNRAS, 349, 722 
Read, A.M., & Ponman, T.J., 2003, A&A, 409, 395 
Sanderson, A.J.R., Ponman, T.J., Finoguenov, A., Lloyd-Davies, E.J., 

& Markevitch, M., 2003, MNRAS, 340, 949 
Sharma, M., McNamara, B.R., Nulsen, P.E.J., et al. , 2004, ApJ, in 

press, astro-ph/0405336 
Sun, M., Forman, W., Vikhlinin, A., et al. , 2003, ApJ, 598, 250 
Sun, M., Forman, W., Vikhlinin, A., et al. , 2004, ApJ, in press. 



|astro- ph/0405511 
Suto, Y., Saski, S., & Makino, N., 1998, ApJ, 509, 544 
Tozzi, P., & Norman, C, 2001, ApJ, 546, 63 
Valageas, P, & Silk, J. 1999, A&A 347,1 
Valageas, P, Schaeffer, R. & Silk, J. 2003, MNRAS, 344, 53 
Valdamini, R., 2003, MNRAS, 339, 1117 
Vikhlinin, A., Forman, W., & Jones, C, 1999, ApJ, 525, 47 
Voit, G.M., Bryan, G.L., 2001, Nature, 414, 425 
Voit, G.M., Bryan, G.L., Balogh, M.L., & Bower, R.G. 2002, ApJ, 

576, 601 
Voit, G.M., Balogh, M.L., Bower, R.G., Lacey, C.G., Bryan, G.L., 

2003, ApJ, 593, 272 
Voit, G.M., Ponman, T.J., 2003, ApJ, 594, L75 
Wechsler, R.H., Bullock, J.S., Primack, J.R., Kravtsov, A.V, & Dekel, 

A., 2002, ApJ, 568, 52 
Wu, K.K.S., Fabian, A.C., & Nulsen, PE.J., 1998, MNRAS, 20, L301 
Xue, Y.J., & Wu, X.P, 2003, ApJ, 584, 34 
Yamada, M. & Fujita, Y, 2001, 553, L145 
Zhao, D.H., Jing, YP, Mo, H.J., & Borner, G., 2003, ApJ, 597, L9 



